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Abstract 

The perturbative quark-mass corrections to the r hadronic width are 
analysed to O(a^rriq), using the presently available theoretical information. 
The behaviour of the perturbative series is investigated in order to assess 
the associated uncertainties. The implications for the determination of the 
strange quark mass from r decay data are discussed. 



1. Introduction 



The inclusive character of the total r hadronic width renders possible an 
accurate calculation of the ratio |1]-§J 

_ r [t~ -> Mr hadrons (7)] 
r [r -> z/ r e i/ e (7) J 

using standard field theory methods. The result turns out to be very sensitive to 
the value of ct s (M%). Moreover, the uncertainties in the theoretical calculation 
are quite small and dominated by the perturbative errors. This has been used to 
perform a very precise determination of the QCD coupling at low energies ||. 

Quark masses play a rather minor role in R T . Owing to the tiny values of m u 
and rrid, their associated corrections are very small || (~ —0.1%). The strange 
quark contribution to the total r hadronic width is suppressed by the Cabibbo 
factor IVusI 2 , which puts the induced m s correction also at the per cent level. 
However, if one analyses separately the semi-inclusive decay width of the r into 
Cabibbo-suppressed modes (i.e. final states with an odd number of kaons), the 
relatively large value of m s induces an important effect of a size similar to the 
massless perturbative correction and of opposite sign |§ . The corresponding R Tt s 
prediction is then very sensitive to the strange quark mass and could be used to 
extract information on this important, and nowadays controversial, parameter. 
A very preliminary value of m s , extracted from the ALEPH r decay data, has 
been already presented in recent workshops ||[7|. 

The determinations of light quark masses are usually obtained from analy- 
ses of the divergences of the vector and axial-vector current two-point function 
correlators or related observables |]8Hl3|. These correlators are proportional to 



quark masses and, therefore, are very sensitive to their numerical values. Unfor- 
tunately, one needs phenomenological information on the associated scalar and 
pseudo-scalar spectral functions, which are not well known at present. The obvi- 
ous advantage of a possible determination of m s analysing quark mass effects in 
r decays is that the experimental error can be systematically reduced in foreseen 
facilities like tau-charm or B factories. There is then some hope to achieve a 
precise determination of m s from such analyses. 

Recently the 0(a^) corrections to the J = quark correlators have been 
calculated |TJ[], and have been found to be rather large. The influence of these 



0(al) corrections on the determination of quark masses and the uncertainties 
coming from the truncation of the QCD perturbative series depend very much on 
the observable. One can see for instance that the QCD perturbative series behaves 
geometrically to 0{a^) for the divergence of pseudo-scalar (scalar) currents if 
resummed perturbatively in terms of a s (s) [T^,[T3|]. This convergence improves 



H| using other resummations like the Principle of Minimal Sensitivity (PMS) |T5| 



or the one advocated in Ref. 
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The hadronic r decay width has also a J = contribution, which, as we shall 
see, behaves rather badly. However, the largest quark-mass correction originates 
in a piece of the left-handed current correlation function, involving the J = + 1 
combination, which shows a much better perturbative convergence. 

The purpose of this paper is to study the perturbative behaviour of the correc- 
tions to R T which are proportional to m 2 ., in order to assess the associated uncer- 
tainties. These are the leading theoretical uncertainties in the m s determination. 
The 0(a s m 2 ) contributions were already studied in Ref. ||. In Ref. |l6j the con- 
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tributions of 0(a 2 m 2 ) to the relevant correlators were worked out. More recently, 
some partial information on 0{a^ s m 2 q ) corrections has become available |14| , |l7| . A 
much detailed analysis of all contributions up to dimension four will be presented 
elsewhere. 

2. Theoretical Framework 

The theoretical analysis of R T involves the two-point correlation functions for 
the vector V£ = ipj^ipi and axial-vector = ipj^^^ipi colour-singlet quark 
currents (i, j = u, d, s): 



n^y(g) = i J rf 4 xe^(0|T(^(x)^(0)t)|0), (2) 
n^ A (g) = • / d*xe^(0\T(A»(x)A^(OnO). (3) 
They have the Lorentz decompositions 

n&v/Ab) = (-<? M V + <f<n n^v/Ai 2 ) + <t<? <,V M (g 2 ), (4) 

where the superscript (J) in the transverse and longitudinal components denotes 
the corresponding angular momentum J = 1 (T) and J = (L) in the hadronic 
rest frame. 

The imaginary parts of the two-point functions n[^ v ^ A (q 2 ) are proportional 
to the spectral functions for hadrons with the corresponding quantum numbers. 
The semi-hadronic decay rate of the r can be written as an integral of these 
spectral functions over the invariant mass s of the final-state hadrons: 



R-r 



127T f**L (l-°) 2 \ll + 2 A;) ImllW (s) + Imll( ) (s) 
Jo Ml \ Ml LI Ml ) 1 ' 1 ' 



(5) 



M 2 ) 

The appropriate combinations of correlators are 

n (J) (.) = \v ud \ 2 (ngvOO + n% A ( s )) + |K S | 2 (ngvW + 0*)) • (6) 

We can decompose the predictions for R T into contributions associated with 
specific quark currents: 

R T = R T ,v + Rt,a + R T ,s ■ (7) 
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R T y and R t ,a correspond to the contributions from the first two terms in Eq. (0), 
while i? Tj 5 contains the remaining Cabibbo-suppressed contributions. 

Exploiting the analytic properties of the correlators n^(s), Eq. (||) can be 
expressed as a contour integral in the complex s plane running counter-clockwise 
around the circle \s\ — M 2 : 



Rr = -irii - 1-4. 3 1 + -^ \D L+T (s)+AD 





We have used integration by parts to rewrite R T in terms of the logarithmic 
derivative of the relevant correlators, 



D L+T (S) EE -8 



d 

ds 



n (0+1) (*)] , DL ^^ipi [ sn(0)(s) l ' (9) 



which satisfy homogeneous renormalization group equations. 

Using the Operator Product Expansion to organise the perturbative and non- 
perturbative contributions to the correlators into a systematic expansion 



in 



powers of l/s, the total ratio R T can be expressed as an expansion in powers of 
1/M%, with coefficients that depend only logarithmically on M T [BJ: 

R T = s(\V ud \ 2 + \V us \ 2 )S EW (l + 8' EW + 5 i0] + E (cos 2 Mi? + sin 2 Mif 

^ D=2A.... 



(10) 

where 5^ = {5\fy + d\f\)/2 is the average of the vector and axial-vector cor- 



rections of dimension D, Sew an d 5' EW contain the known [T9],[2(J electroweak 
corrections, and sin 2 ^ ee |K s | 2 /(|Kd| 2 + \V US \ 2 ). 

The dimension-zero contribution is the purely perturbative correction neglect- 
ing quark masses, which, owing to chiral symmetry, is identical for the vector and 
axial-vector correlators. It is fully generated by the Adler function D L+T (s), be- 
cause D L (s) vanishes in the chiral limit. The correction 8^ has been investigated 
in great detail in Ref. [Q. We will follow a similar procedure to analyse the per- 
turbative quark-mass corrections of dimension two. 

3. Dimension— Two Corrections 

For the sake of simplicity, let us take here m u = = 0. The arguments we 
shall put forward don't depend on it. In this limit, the vector and axial- vector 
correlators get the same quark-mass corrections, i.e. D^ s A (s) = D^ sV (s) = 
Du S (s) (J = L + T, L). The dimension-two contributions can be written in the 
form: 

D L u f{s)\ D = 2 = ^ £ d L n +T (0 a n (-es) , (11) 

'* 7r 6 n=0 

D L us (s)\ D=2 = - A^l^ E^(0a"(-6), (12) 

° 7 ' 1V1 T 71=0 
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where a = a s /ir, £ is an arbitrary scale factor (of order unity) and the coefficients 
d-niO are constrained by the homogeneous renormalization group equations sat- 
isfied by the corresponding functions D^ s (s): 



£ 4 dt(0 = E [2 7fc - (n - k)0 k ] d J n _ k (0 

' k=l 



(13) 



for n > 1 and 



i.e. 



(14) 



#(0 



d J 



< + 2 7l ^log£. 



di + 



2 72 < + (2 7l - Pi)d{ log£ + 7l (2 7l - (3M log 2 £ . 



di+ 2 73 ^ + (2 72 - 2 )d( + 2(7! - log£ 



(15) 



+ 



(- 7l /3 2 + 2 72 (2 7l - /3i)) d J + ( 7l - A)(2 7l - ft)<tf log 2 £ 



+ 3 7 i( 7 i-A)(2 7l -AK J log 3 £, 

<S(0 = <# + ••• 

The factors and 7 fc are the expansion coefficients of the QCD and 7 functions, 



= /5(a) a. 

d/i 



(3(a) = 



(16) 



fc=i 



dm 
d/i 



-7(a) m. 



7(a) = E7fc flA 



;i7) 



/c=i 



which are known to four loops |2l|-p3l. The coefficients d J n = d^(l) are only 
known to order a 2 for J = L + T _and a 3 s for J = L [§,[g,[i^|2|-[28|. 
For three flavours and in the MS scheme, one has: 



0i 
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04 



02 = 

140599 445 



2304 



16 



71 



2; 

7 4 



72 =12 



2977517 
20736 



91 
12 
9295 
216 



03 = - 

-94.456 079 
8885 



3863 
" 192 ' 



73 



. 135 
Ca + -=r Ca 



288 
125 
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C5 



5Cs « 24.840 410 ; 
! 88.525 817 ; (19) 
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iL+T 



1 ; 



jL+T _ 13 . 

1 3 ' 



lT , T 21541 323 _ 520 „ 

4 = + Cs =- C5 ~ 37.083 047 ; 

2 432 54 s,j 27 s 



, T 9631 35 _ 

do = Cs 

2 144 2 ss 



45.845 949 



(20) 



, r 4748953 91519 5 , 715 , 

4 = C3 C4 H Cs 

3 5184 216 s 2 s 12 s 



465.846 304. 



Notice the rather bad perturbative behaviour of the D = 2 corrections to the 
correlation functions D^ s . Remember that a(M 2 ) ~ 0.11. 

Inserting the expansions (|ll]) and (|12|) in Eq. (j8|), the D = 2 corrections to 
R T> s can be expressed as 



6 (2) 
w us 

where 



Ml 



A[a(M T 2 )] 



AH^lfsA^H + A^a]}, (21) 



A J [a(M 2 )] = E^)4 n) W, 



(22) 



n=0 



and the contour integrations are contained in the functions 



-1 



47TZ J Lcl=l 3^ 

1 /" cfe 

27TZ 7|a|=l X 



X) 



m(M T 2 ) 



a"(-£ 2 M 2 x), (23) 



m(-^ 2 M 2 x 



(24) 



Since the quark mass ratio is flavour independent, the integrals B < f\a^) regulate 
also the small corrections proportional to m u and m^, which we are neglecting. 
These functions depend only on 7j -, = a s (£ 2 M 2 )/n and log£. Moreover, 
they satisfy the following homogeneous renormalization group equations: 



(25) 



k=l 



4. Perturbative Expansion 



The usual perturbative approach expands the Bj(a^) functions in powers of 
d£. This gives, 



J Bi 0) (a c ) = l- 7 i 2\ogi + H( 



272 io g e - 71(271 + a) ^g 2 e + (72 - 2 7l 2 iog£) 
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2 73 log£ - (2 72 (2 7l + 7l /3 2 ) log 2 £ + -71(71 + A) (27i + A) log 3 ^ 



+ (73 - 4 7 i72 logC + 7i 2 (27i + A) log 2 H( 

+ ^i72 + |&7i - 7i72 + 7? [li - y) log^ #1 



+ ^(2 7l -/? 1 )(7i-A)^ 3 J 



a| + 



(26) 



5? ) (a 5 ) = a 5 - 



2 7l log£+f 7l -^W J 



^2 



2 72 log(0 - 71(271 + A) log 2 £ + ( 72 - y - 7i (27i - A) loge) H( 



a\ + 



4 2) (a 5 ) = a\ - |2 7l log£ + (71 - A) 



a| + 



B< 3) (a € ) = af 
where 



Hi +T ^jii. .^(1 + »)(!-* ) 3 log" (-x), 



— 1 /" 

4m /|x|=i x 2 

^^ii-if (1_l)3log " (_l) - 



To 0(etg), the needed integrals are 



H L+T = 


1, 






1 

~ 6 ' 


Ht 


_rp 25 

~~ 18 


7T 2 

3 




85 
" 36 ~ 


7T 2 

6 


1 


# 4 L+T = 


721 

54 


25 

— IT 2 + 

9 


5 ' 


H L = 


1, 




H? = 


11 




85 
" 18 ~ 


7T 2 

3 ' 


H 3 L = - 


575 
36 + 


11 
~6 


vr 2 , 




3661 
54 


85 2 
9 
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The perturbative expansions A J [a] then take the form 



A J [a(M T 2 )] = J2[dn(0 + hn(0 



n=0 



(27) 
(28) 
(29) 

(30) 
(31) 



(32) 



(33) 



(34) 
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where the coefficients h^(^) depend on c^ <n (£), f3 m<n and 7 m < n ; thus, they are 
known up to 0(a 3 ) and 0(a 4 ) for J = L + T and J = L respectively. For £ = 1, 
one has [ft^ = ^(1)]: 

4 +T = ; h^ T = - 1 - ; h L 2 + T = 1 -§~ 1 ^ 2 ^ -12.412 495 ; 

r , T 114517 4391 , 3659 , 1690 , , . 

hk +T = vr 2 Cs + ~ -241.926 329 ; 35 

3 2592 144 648 S3 81 S5 ' V ; 

lL+T 26864009 3110783 2 1073 4 3051761 . 123745 2a 

ft 4 = 7T 7T H Cs Cs 

4 13824 5184 1920 15552 s 2592 s 

171845 , 29575 2a 35 lT , T mnm 35 lT , T 

Cs + —T^T T Cs - 777 4 ~ ~ 3 229.101 787 - — 4 ; 

243 162 s 24 3 24 3 

Jtf = 0; ^ = Y' ^ = ^-11^^ 64.143 060; 

, T 1435691 7927 9 5225 M , . 

4 = vr 2 Cs ~ 856.673 579 ; 36 

3 864 144 24 S3 ' K J 

L _ 693706385 2295071 2 10877 4 1429525 5595 2 264275 

4 ~ 20736 1728 ^ ~ 17280 * 144 ^ 3 + ~32~ * ^ 3 + 288 ^ 5 

« 11377.111254 . 

The contour integration generates rather large numerical factors, which show 
an opposite behaviour for the transverse and longitudinal pieces. In the A L+T ex- 
pansion the h^ +T contributions cancel to some extent with the original correlation- 
function coefficients d^ +T , 

A^[a(M 2 )] = 1 + 4 a(M 2 ) + - ^ + ^ C s - ™ Cs) «(M 2 ) 2 

/,, +T 114517 4391 2 3659 , 1690 \ 

= 1 + 4 a(M 2 ) + 24.671 a(M 2 ) 2 + (4 +T - 241.926) a(M 2 ) 3 + • • • (37) 

However, both d% and /i^ contributions are large and positive, which gives rise 
to a badly behaved expansion for A L : 



/ 13363099 7927 , it 4 17318, 715 , \ ,,,,,, 
+ (-6184 - " UA * ~ 36 " — 6 + U Cs J a(M; » + " ' 
= 1 + 9.333 a(M 2 ) + 109.989 a(M 2 ) 2 + 1322.520 a(M 2 ) 3 + • • • (38) 
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The bold-guess estimate d\ ~ d\ (d^/d^) ~ 4733 would result in a huge 0(a 4 ) 
coefficient d\ + h\ = 16 110. 

Since A L+T has a larger weight on the total contribution to 5^], the final 
combination of the transverse and longitudinal pieces has a better behaviour^: 

W)] = 1 + f a(MS + (i5™ - II ^ + i (3 - f ft) <Mif 

(3 14050201 _ 5275 _ ^ _ 47401 13195 \ 

\4 3 20736 144 144 288 432 s ) y r> 

+ ••■ 

= 1 + 5.333 a(M 2 ) + 46.000 a(M 2 ) 2 + fl49.185 + ^4 +T ) «(M 2 ) 3 + • • • 

(39) 

Nevertheless, the convergence of this perturbative series is very poor for the 
range of the strong coupling relevant in r decays, a(M 2 ) ~ 0.11. With d^ +T ~ 
g?2 +T (^2 +T /^i +T ) ~ 317, the 0(a 3 ) correction would be of the same size as the 
0(a) and 0(a 2 ) contributions. 

5. Resummation of Running Effects along the Integration Contour 

At the moment, we can do very little about the apparent growth of the d 3 n 
coefficients, specially for J = L. We clearly need a deeper understanding of the 
perturbative D 3 s (s)\ D=2 expansions. However, we can try to control better the 
large contributions contained in the h 3 n factors. 

The integration along the circle x = gives rise to a long running of the 
quark mass and the QCD coupling. The expansion of m 2 (— £ 2 M 2 x) a n (— £ 2 M 2 a;) 
in powers of generates imaginary logarithms log™ (— x) = i n (<j) — 7r) n , which are 
large in some parts of the integration range. The radius of convergence of such 
expansion is actually quite small Q. However, there is no need to perform this 
ill-defined power expansion. 

Using in Eqs. (|23"D and ( p4[) the exact solution for m(— s) and a(— s) obtained 
from the renormalization group equations, the Bj l \a^) integrals can be calculated 
to all orders in a s , apart from the unknown /3 n> 4 and 7 n >4 contributions, which 
are likely to be small. Thus, a more appropriate approach is to directly use the 
expansions (|22|) , in terms of the original d 3 n coefficients, and to fully keep the 
known four-loop information on the functions Bj l \a^). 

Tables [1] and show the exact results for i?j"^ T (a) and Bj^\a) (n = 0, 1, 2, 3) 
with £ = 1 obtained at different orders in the (3 and 7 expansions, together with 



a The 0(a 2 ) correction agrees with the numerical result recently reported in Ref. |29(] , which 
is larger than the value originally quoted in Ref. fill . This larger 0(a 2 ) correction has been 
also confirmed by K. Chetyrkin and A. Kwiatkowski pQ). 
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the final values of A J [a], for a = 0.1 (£ = 1). For comparison the numbers coming 
from the truncated perturbative expressions at 0(a 3 ) are also given. 

Table 1: Exact results for B^_ T (a) (n = 0, 1,2,3) obtained at the k-loop (k = 
1,2,3,4) approximation (/3j>k — lj>k — 0), together with the final value of 

A L+T [a] = En=o rf n +TjB i+T( a )> for a = 0.1 and £ = 1. For comparison the 
numbers coming from the truncated expressions at 0(a 3 ) are also given. 



Loops 








B^ T (a) 


A L+T [a] 


1 


0.890 32 


0.069 65 


0.004 52 


0.000186 


1.360 


2 


0.81719 


0.056 66 


0.002 78 


-0.000 008 


1.166 


3 


0.79143 


0.052 96 


0.002 36 


-0.000 048 


1.108 


4 


0.782 37 


0.05168 


0.002 22 


-0.000 060 


1.089 


0(a 3 ) 


0.793 63 


0.064 73 


0.008 92 


0.001000 


1.405 



Table 2: Exact results for (a) (n = 0,1,2,3) obtained at the k-loop (k = 
1,2,3,4) approximation (/3j>k — lj>k = 0), together with the final value of 
A L [a] = Y^ n =Q d^B^\a) , for a = 0.1 and £ = 1. For comparison the numbers 
coming from the truncated expressions at 0(a 3 ) are also given. 



Loops 


Bf{a) B [ l\a) Bf\a) B { ^{a) 


A L [a] 


1 
2 
3 
4 


1.399 08 0.184 73 0.022 55 0.002 588 
1.54013 0.202 47 0.024 21 0.002 692 
1.578 53 0.20617 0.02444 0.002 690 
1.58910 0.20706 0.02446 0.002 681 


4.686 
5.052 
5.120 
5.133 


0(a 3 ) 


1.64446 0.218 94 0.02192 0.001000 


4.356 



These numerical results show a reasonable convergence of the Bj (a) inte- 
grals, as higher-order /3k and 7^ contributions are taken into account. Increasing 
the number of loops one gets a small decrease (increase) of the transverse (lon- 
gitudinal) contribution. It is also clear that the truncated 0(a 3 ) expressions 
overestimate (underestimate) A L+T (A L ). Taking the full four-loop information 
into account, we get the following perturbative behaviour: 

A i+T [0.1] = 0.7824 + 0.2239 + 0.0823 - 0.0000601 d^ +T + ■■■ (40) 

A L [0.1] = 1.5891 + 1.1733 + 1.1214 + 1.2489 + • • • (41) 

The L+T series converges very well. Owing to the negative running contributions 
the A i+T [a] series behaves better than the original perturbative expansion of 
^ > us~ T ( s )\d=2- Unfortunately, the longitudinal series is much more problematic. 
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The bad perturbative behaviour of D^ s (s)\ D=2 gets reinforced by the running 
effects, giving rise to a badly defined series. 
The combined final expansion, 

A[0.1] = 0.9840 + 0.4613 + 0.3421 + (0.3122 - 0.000045 4 +T ) H (42) 

looks acceptable for the firsts terms because 

A L+T 

is weighted by a larger factor. 
In fact, this series behaves better than the one in Eq. (|3"9"|), obtained with the 
usual perturbative truncation of the contour integrals. Nevertheless, after the 
third term the series appears to be dominated by the longitudinal contribution, 
and the bad perturbative behaviour becomes again manifest. 

Using the full four-loop result we have certainly gained in convergence for the 

A L+T 

series [compare the fourth term in the series ([J7]) and ( |40"D for a = 0.1], 
which is otherwise the one we don't know the 0(a 3 ) coefficient. We can take 
advantage that the 0(a 3 ) correction to A [a] is almost completely given by the 
known A L contribution. Using d^ +T ~ d 2 +T (d 2 +T '/ 'd^ +T ) « 317, the fourth term 
in P^) becomes 0.298, i.e. a 5% reduction only. Taking the size of the 0(a 3 ) 
contribution to A L as an educated estimate of the perturbative uncertainty, we 
finally get 

A[0.1] = 2.1 ±0.3 . (43) 

6. Renormalization— Scale Dependence 

The expansion fl22| ) depends order by order on £ and this dependence cancels 
out only when we sum the infinite series. In practice, we only know a few first 
terms of the series (three for A L+T [a] and four for A L [a]); so we should worry 
how much the predictions depend on our previous choice £ = 1. Obviously, £ 
should be close to one in order to avoid large logarithms; but variations within a 
reasonable range, let us say from 0.75 to 2, should not affect too much the final 
results. Smaller values of £ would put the QCD coupling in the non-perturbative 
regime and are therefore not acceptable. 

Figures [I], |2| and |3| show the sensitivity to the selection of renormalization 
scale of the final predictions for A L+T , A L , and A, respectively, for a(M^) = 0.1. 

The behaviour of A L+T is quite good. The predicted value remains very stable 
in the whole range £ G [0.75,2], showing that the perturbative series is very 
reliable. Below £ ~ 1/2, the perturbative expansion breaks down, as expected, 
because the coupling is already outside the radius of convergence of the series. 

The longitudinal series, on the other side, has a quite wild dependence on the 
renormalization scale. Changing £ from 1 to 2, amounts to a reduction of A L of 
about 65%. Thus, the theoretical uncertainty is very large in this case. 

The £ dependence of the complete expansion A, reflects obviously the be- 
haviour of its two components. The larger weight of A i+T keeps the result still 
acceptable, within the range of £ considered, but the sizeable A L contribution 
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spoils the stability and generates a monotonic decrease of the prediction for in- 
creasing values of £. Taking this variation into account, the theoretical error in 
Eq. (|43|) should be increased to about 0.6, i.e. a 30% uncertainty in the final 
prediction. 



6. Discussion 



The bad perturbative behaviour of the longitudinal contribution does not 
allow to make an accurate determination of the strange quark mass from R Tt s- 
Nevertheless, taking 

A[0.1] = 2.1 ±0.6 , (44) 

m s (M^) could be still obtained with a theoretical uncertainty of about 15%, 
which is not so bad. 

Notice that it is the phase-space integration of the original correlation func- 
tions the responsible for the different behaviour of the longitudinal and trans- 
verse components. Therefore, the perturbative convergence could probably be 
improved through an appropriate use of weight factors in Eqs. (§) and (||). This 
requires an accurate measurement of the final hadrons mass distribution in the r 
decay, which so far has only been performed for the dominant Cabibbo-allowed 
modes [31]. The measurement of R us (s) could be feasible at the forthcoming 
flavour factories, where a very good kaon identification is foreseen. 

From the theoretical point of view, the analysis of weighted moments of the 
final hadrons mass distribution proceeds in a completely analogous way A 
detailed study will be presented in a forthcoming publication. 
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Figure Captions 



• Figure 1.- Variation of Ai + ^[0.1] with the renormalization-scale factor £, 
to four loops. 

• Figure 2.- Variation of A L [0.1] with the renormalization-scale factor £, to 
four loops. 

• Figure 3.- Variation of A [0.1] with the renormalization-scale factor £, to 
four loops. 
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Figure 1: Variation of A i+T [0.1] with the renormalization-scale factor £, to four 
loops. 



A L [0.i] 

10 i 1 1 , 




Figure 2: Variation of A^[0.1] with the renormalization-scale factor £, to four 
loops. 
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Figure 3: Variation of A[0.1] with the renormalization-scale factor £, to four 
loops. 



